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POSSIBILITY OF A NEW APPROACH TO PROCESSING SPACE 
PULSE RADIO SIGNALS 

A. A. Bocharov 


Proposed herein is a new approach to the processing /2* 

of space pulse radio signals, which makes it possible to 
achieve a high level of time resolution. The approach is 
based on the utilization of the "Fourier transform of a 
radio pulse envelope in the interstellar mediurA" effect. 

The conditions of applicability of the proposed method 
and the requirements for recording equipment are substan- 
tiated. Also i^:?,dicated is the possibility of obtaining 
"superresolution'' with the utilization of the givjn 
method of signal processing. 


Introduction 


During studies of the internal structure of pulse signals of 
cosmic origin, for example, pulsar radiation, substantial limitations 
are placed on time resolution by effects associated with the propa- 
gation of radio waves in the interstellar plasma,— .with the dispersion 
and scattering on discontinuities. The dispersed state of the 
interstellar medium leads specifically to distortion of form and 
broadening of the pulses, in connection with which, in order to 
ensure a high level of time resolution, the development of special 
methods of "elimination of dispersion" [3] was required, which lead 
to an appreciable complication of both the recording equipment and 
the subsequent digital processing of the received signals. 

However, as will be shown below, during the implementation of 
certain relationships between the parameters of the radiated 
pulses and the interstellar medium, another approach is possible to 
processing of the space pulse radio signals, with which the dispersed 
state of the medium may prove to be a "useful" effect, which makes 
it possible to study the internal structure of the pulses with a 
high level of time resolution, with a simultaneous reduction in 
the requirements for the recording equipment. The proposed method, 

♦Numbers in the margin indicate pagination in the foreign text. 
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which does not require comp, lex specialiised equipment and lengthy 
digital proces'iing of the received signals, may, in our opinion, 
find application in studies of pulsar radiation, and possibly for /4 
problems of space communications as v^ell. 

1 . Connection Between Emitted and Received Signals 

We will mention at once that the entire subsequent examination 
will be* carried out on the assumption of the narrow-bandedness of 
the examii>ad signals, which is fully justified, insofar as in 
radioastronomy, as a rule, the band width of the receiving apparatun 
is small, as compared with the central frequency. The emitted 
signal, like a narrow- band signal, may be written in the form 

« K (*) I 

and, accordingly, the received signal may be written as 

iiM - C0i> (U^ot (t)) , 

where Ai ^(t) and 2^t) are the envelope and phase of the signals, 
wo is the carrier frequency (for broad-band pulsar radiation, wq 
is understood to mean the central frequency of the receiver) . 

In order to determine the connection between Si(t) and S2(t), 
it is advisable, in view of their narrow-bandedness, to utilize 
the method of complex envelopes [ 1 ], Let ^i(t) and ^he 

complex envelopes of the signals Si (t) and S2(t), determined by the 
expression 

The emitted and received signals are associated by the xelationship 

^ y ( 1 ) 

where Gi 2(w) are the spectra of the complex envelopes di 2 (t) f 

, f ' 

H(w) is the equivalent low-frequency coefficient of transmission of 
the interstellar medium. 


2 


For narrow-band signals, in the approximation of geometric 
optics ( i»e. , with the condition that tht parameters of the medium 


change little on the wave length, which i« indisputably correct 
for the frequencies utilized in radioastronomy), i)(u) has the 
form [3], [51 



( 2 ) 


(In expression (2) for , two terms are omitted, which are 
insignificant for subsequent examination — the constant phase 
shift, which is lost during detection of the envelope, and also 
the term which determines the time lag between the moments of 
radiation and receiving of the pulse). The parameter olq in (2), 
which has the meaning of the derivative instantaneous frequency 
of the received signal according to time, is determined, in turn, 
as 


d 


0 
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where C, eq, m, e are the speed of light, the dielectric permeability 
of the vacuum, the mass and charge of the electron, respectively; 

<iv ^ 

^/M=o/ Ng(z)dz is the degree of dispersion of the radio source, 

which is the integral of the electron concentration along the 
trajectory of propagation of the radio wave. 


According to (1) and (2) , the expression for the complex 
envelope of the received signal takes on the form 
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2. stationary Phase Method and Conditions of its Applicability 
The stationary phase method (MSP) is widely utilized for 
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approximate calculation of integrals of the type (3) [4]. Assume 
that we have the integral 

I = J </<*/ , 

where 4>(u) is the actual function, and R((») may also be complex. The 
stationary phase method is applicable in that case when ♦(u) contains 
even more rapidly changing components of the phase, and R(ai) is 
a slowly changing function. This shows that the basic contribution 
to the integral, in this case, is provided by the points of 
"stationatiness" of the function 4>(aj), in connection with which 

IS R ran V [• ^ exp [- ) 

where u)g is the point at which 4>(uj) is stationary, i.e. , =0. 

S 

With the presence in the function 4>(a)) of several stationary points, 
it is necessary to sum up the contribution of each of them. 

For simplicity, we will take a radio pulse without phase 
modulation as the emitted signal, i.e. , we assume that Vi(t)=0. 

(The examination with regard for the presence of phase modulation 
in the signal will be carried out below) . We will rewrite (3) 
in the form 

%(*) “ zr) K e ‘ 

where Wi (w) is the spectrum of the envelope A_(t) of the emitted 

2 ^ 
radio pulse, and F(a))=“( 5 — +wt) . 

The problem consists of finding the condition of applicability 
of the stationary phase method for calculating the integral 
in (4). Comparison of the functions Wi (w) and will be 

carried out with t==0, which corresponds to the case of the "slowest 
change" in the function F(w,t). Let the chara<?teristic width of 


the emitted pulse be equal to then, accordingly, the 

characteristic width of its spectrum (which contains the basic 
portion of the pulse unergy) '^“en“^’^'^^pulse* order for the 
function Wi((o) to be assumed to be slowly changing# as compared 
to , it is natural to adopt the following condition 

(see Fig. 1) : 



Let there be a large number of "periods” of the function 
in (here# "period" is understood to mean the distance between 

©n — ip/fii) 

two adjacent points# in which the function e ' ' takes on equal 
values# for example, equal to I, which takes place with F(o))=2Trn, 
where n is an integer^ . The number of these "periods" in 
will be equal to F(u)=Ato )/2ir# and, consequently, the condition of 
the slow change in the function Wi (w) , as compared with e 
written in the form 


-iF (u) 


IS 
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or, with regard for the connection between and Ato 


_L 

ztr 


en 


o(-t* , << 'H' 

'^0 pulse 


In practice, in the majority of cases, one may utilize the 
following condition of applicability of the stationary pulse 
method, with which we will subsequently operate, namely 





( 5 ) 


■<kj\ 4 1 . 

0 pulse 

3. Fourier Transform of the Envelope of Pulse Radio Signals in 
the Interstellar Plasma 


Let the parameters of the pulse and the medium satisfy condition 
(5). By calculating the Interval In (4), using the stationary 
phase method, and taking into account the fact that, in the given 
case, the function F(w) has only one stationary point Ug*-aot, 
we obtain the following for the complex envelope of the received 
signal 



Wth***-) ‘ 6 ^ 



The "physical” envelope of the received signal ( i.e. , the voltage 
at the output of the linear detector) will hava the form 


, 6 , 

(It is taken into account here that Wi (w) is the spectrum of the 
actual function Ai (t) , and, consequently, jWi (-to) | = |Wi (w) | ) , 


Ignoring the change in scales, one may establish that, with 
the implementation of condition (5) , the envelope of the received 
signal is the modulus of the spectrum of the envelope of the 
emitted signal, i.e. , the interstellar medium has accomplished its 
own type of Fourier transform of the envelope of the initial radio 
pulse. 

By squaring both parts of (6), we will obtain the following for 
the square of the envelope of the received signal ( i.e. , the 
voltage at the output of the quadratic detector) 

( 7 ) 
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The magnitude o£ |Wi(ci))p Is, in essence, the energy spectrum 
o£ the signal Ai(t). We will find the inverse Fourier transform 
from both parts of (7) , t 2 Ucing into account the f 2 ^ct that, 
according to Fourier, the autocorrelation function (ACF) of the 
signal corresponds to its energy spectrum. We have 

Xit) * ^ K(^) (8) 

+ 00 

where / Ai (t) Ai (t-fT)dt is the autocorrelation function of 

the envelope Ai (t) , and X(t)— 3C(t)e^^Mt. 

p-00 * 

Thus, with the observance of condition (5), the inverse 
Fourier transform from the square of the envelope of the received 
signal will give us the autocorrelation function of the envelope 
of the emitted signal. What is more, as will be shown below, the 
presence in (8) of the multiplier I/ocq in the argument of the 
autocorrelation function makes it possible to obtain the auto- 
correlation function of the envelope of the initial pulse with 
"superresolution", i.e . , with a time resolution higher than at 
the recordi«g equipment. 

4, Basic Time Relationships 

As was shown above, with the observance of condition (5), for 
obtaining the autocorrelation function of the envelope of the 
emitted radio pulse, it is sufficient to find the inverse Fourier 
transform from the square of the envelope of the received signal. 

We will now determine the time resolution obtained utilizing the 
given method of processing of the signals, and the requirements for 
the recording equipment. (Everywhere below, where it is not 
especially specified, we will be implying requirements for the 
low-frequency — i.e. , after the detector of the envelope or 
heterodyning to low frequencies — portion of the equipment. In 
practice, with the utilization specifically of the method of 
"elimination of dispersion", it is namely the low-frequency portion 
of the receiving-recording equipment. Which includes a rapid-action 
analog-digital converter, external information storage units with 
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a device for controlling them, and so on, which is the fcaaic 
complexity and determines, in the final analysis, the time 
resolution of the entire system) . 

Let T be the observation interval, At is the quantization 
spacing (time resolution) of the recording equipment, and N«T/At 
is the number of readings of the envelope A 2 (t) of the received 
signal, fixed during the interval T, After the inverse Fourier 
transform of the square of the envelope, we obtain the signal X(t), 
determined in the interval O' with a quantization spacing At', 
with 


According to (8), the quantization spacing At' of the function 
X(t) corresponds to the fact that we obtain the values of the 
autocorrelation function of the envelope Ax (t) with a quantization 
spacing AT=AT'/ao. We also take into account the fact that the 
discrete inverse Fourier transform of the function X(t) will givs 
us the function X(t), which is symmetrical relative to the middle 
of the interval 0^ , i.e. , in actuality, we have only N/24-I 
independent readings of the autocorrelation function [2], Thus, 
we obtain the autocorrelation function of the envelope of the 
initial radio pulse with a quantization spacing 


(9) 



and determined in the interval 




At 


In order to completely obtain the autocorrelation function 
(and,, what is more, avoid undesirable effects associated with 
cyclic convolution) , it is necessary to fulfill the condition 
0>tpuige/ from which we. obtain the requirement for the time 
resolution of the recording equipment 

At i j 


’‘•h 


ulse 


( 10 ) 
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The number of readings of the envelope of the rec,;sived signal, 
which must be fixed during the interval T, is accordingly equal 
to 



f 


We will now explain the meaning of the term ’•superresolution”, 
which has already been utilized. Let the equipment resolution 
be selected as the “worst'' permissible accn;rding to (10) , i.e. , 

At = -2^ 

*:;?hen, from (9) and (11), we obtain the equation of "super- 
resolution" 


Consequently, taking the observation interval as sufficiently 
large, as compared with the duration of the emitted pulses, we 
may obtain the autocorr^ilation function of the envelope Aj (t) with 
a time resolution considerably better than in the recording 
equipment. (In practice, the selection of a value of T which is 
as large as one pleases, and, consequently, obtaining a resolution 
as high as one pleases, is impossible because of the presence of 
noises) . We would also note that the band of the receiver (prior 
to the envelope detector) should be matched with the magnitude 

t thus, the term "superresolution" has 

meaning only as applied to the low-frequency portion of the 
recording equipment. 


/II 


We will now evaluate the best resolution, which may be obtained 
utilizing the given method, for pulsars as periodic radio sources. 

It is evident that, for this purpose, one should (with the absence 
of noise) select an observation interval equal to the period of 
the pulsar P, i.e. , The best resolution, in that case, 

will be 


A,P U2) 
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Taking according to (11) , w© obtained a nuirdber of readings 

of the envelope hz (t) , which must be recorded during the interval 
T 



7T 


(13) 


The equation of "superresolution” will now have the form 

■ A V = 2. 

» 

Insofar as, for all pulsars, tpuj^gg/P« 1, consequently, ^topt* 


For each space pulse radio source at a fixed frequency 

£q^(hd/2k, one may introduce the parameter max' which is the 

maximum duration of the emitted pulse, with which condition (5) 

is still fulfilled, i.e., t^. , *l/an. The values of t , 

r pulse max voiuco uj. i-puj^gg 

for some values of 3)h and fo are given in the Table. 


MAXIMUM PULSE DURATION 


1 SM [ps/ 


Is 

• 

iio 

• 

1 20 

• 

■■■■wawiK >txA 

i 50 

V 

j 100 
• 

1 200 
• 

j 500 

^ pulse 

max 

Tinsoc J 








( fu = 

100 

NUIa ) 

0«08 

0,11 

0,16 

0. Z 5 

0,35 

0,5 

0,8 

pulse 

max 

C msec 3 








( To * 

40 

MIJa ) 

0,32 

0,45 

0,64 

I.o 

1,4 

2.0 

3,2 

t>pulse 

max 

r msec ] 








( 

20 

MUz ) 

0,9 

1.3 

1,8 

2,8 

4.0 

5,6 

9,0 


Equation (12) may then be rewritten in the form 
and relationships (11) and (13) will be written as 
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( 14 ) 


4 -ST — “ 

” ''pulse max 

with the eign in (14) occurring with tp„x,e“*pulse mex 
corresponding to the "worst” case, i»e. > the highest requirements 
for the recording equipment. 

5, Comparison with the Method of "Elimination of dispersion" 

We will now examine some advantages and shortcomings of the 
proposed method for processing pulse radio signals of cosmic 
origin. We will carry out the comparison with the aforementioned 
method of "elimination of dispersion" of Hankins (which is, 
in essence, inverse filtration of the received signals, i.e . , 
multiplication of their spectra by the function which is inverse 
to the coefficient of transmission of the interstellar medium) , 
which makes it possible to achieve considerably higher resolution 
than other processing methods. 

The two basic shortcomings of the proposed method are evident 
— first, the impossibility of recreating the very form of the 
envelope of the initial radio pulse (we obtain only its autocorre- 
lation function) ; second, an entire series of restrictions on the 
applicability of the given method, which may be utilized only 
with fully determined relationships between the parameters of the 
pulse and the medium. The Hankins method is completely devoid of 
both these shortcomings. 

The basic advantage of the processing method examined above 
Consists of the considerable reduction in the requirements for the 
recording equipment (more precisely, as was noted above, for its 
low-frequency portion) . Let At be the obtained resolution of the 
autocorrelation function of the envelope of the initial signal 
(or the very envelope for the Hankins method) , and At is the time 
resolution of the equipment. With utilization of the procedure 
of "elimination of dispersion", the obtained resolution of the 


envelope, in the best case, is equal co tbat of the equipment, 
i.e, , AT>At, whereas, for the proposed processing method 
according to the equation of "superresolution”, the requirements 
for the equipment resolution are reduced by times, 

which, in practice, may comprise two or more orders. 

The advantages of the given method may include the fact that, 
in this case, it is not necessary to know the degree of dispersion 
of the radio source with very high accuracy, insofar as it 
determines only the time scale of the attained autocorrelation 
function of the envelope; at the same time, utilization of the 
procedure of "elimination of dispersion" requires the knowledge of 
the degree of dispersion with very low absolute error. /I 4 


It makes sense to dwell here on the matter of the effect of 
scattering of radio waves in the interstellar medium. Scattering, 
like dispersion, leads to broadening of the pulses, which will be 
characterized by the magnitude of '’'^gcatt' scattering 

time constant. We would note that, insofar as there is no 
functional dependence between the degree of dispersion and the 
scattering time constant, the magnitude of '^scatt ® specific 
radio source at a given frequency may only be determined 
experimentally. 


Scattering places a natural restriction ct the time resolution 
of the received signals, which, irrespective of the parameters of 
the receiving equipment, may not be better than the scattering 
time constant. Consequently, with utilization of the method of 
"elimination of dispersion", it is impossible to achieve a time 
resolution better than the magnitude of same 

time, the proposed method, which places requirements T/2tp^j^gg 
times lower on the time resolution of the received signal, makes 
it possible, accordingly, to achieve resolution several times 
greater than the magnitude of In this case, however, the 

scattering places an evident restriction on the applicability of 
the given method, and namely, ■^sc^tt-'^^opt' ^^® roagnitude 

of At^pt is determined from expression (11) , i.e, , 
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r < jai- 

*'scatt ^ ^ + 

**0 *'oiilse 

6. General Case of Radio Pulses with Phase Modulation 


Examined above was the case whei. phase modulation is absent 
in the emitted pulses; we will now make an attempt to evaluate to 
what consequences its presence loads. In this case, for the 
complex envelope of the received signal, we have 


where 


-» i F(^K 

3z(t) ^ ^ ^ 

=* / (*) J F(o>1~ “ (zXl 

^09 


) 


(15) 
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Here, it is difficult to clearly formulate the conditio.i 
of applicability of the stationary phase method for the general 
case without any information on the law of phase change. In the 
extreme case, one may assume that, insofar as the spectra of a 
radio pulse of one and the same form, with the presence in it of 
phase modulation, will be broader than in the case of its absence, 
and, consequently, the characteristic width of the spectrum of 
the complex envelope Gi (u) should be greater then the spectrum of 
the ’’physical'' envelope Wi (co) , the function Gi (w) is "slower 
changing" than Wi (w) , and one may utilize condition (5). But, 
strictly speaking, it is impossible to recognize similar discussions 
as correct, and condition (5) here is inapplicable in the general 
case . 


Let the condition of applicability of the stationary phase 
method be observed. By calculating the integral in (15) using 
the stationary phase method, we obtain the following for the complex 
envelope of the received signal 



and, accordingly, for the square of the "physical" envelope 
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(16) 


x(t ) = ztidt - r 

As is evident from (16), the square of the envelope of the received 
signal is proportional to t >e energy spectrum of the complex 
envelope of the emitted signal, but from a negative argument. The 
energy spectrum and the autocorrelation function of the complex 
signals are associated between themselves by a Fourier transform, 
i.e. , 

X(t) » ir i ( 17 ) 

«• # 

where jGi (a))| ^ is the energy spectrum of the complex signal 0 i (t) ; /16 

and J^(t)=/fl*(t) rtj(t+T)dt is its autocorrelation function. In 

addition, the function J<(-t) will correspond to the function 
1 Gx(-w)P, according to Fourier, with J^(-t) = ^*(t) . 

From (16) and (17) , we obtain the following expression for the 
function X(x), which is the inverse Fourier transform of the 
function.X(t) — the square of the envelope of the received signal: 

X(r) 

XV) = iFxl<(£) . 

In the given case, consequently, we obtain the autocorrelation 
function of the complex envelope of the received signal, and not 
its "physical" envelope, like in the case of pulses without phase 
modulation. One may establish that such a complex autocorrelation 
function for pulses with phase modulation is more informative 
than simply the autocorrelation function of the envelope, insofar 
as it carries information on the law of phase change. The time 
resolution of the obtained autocorrelation function, and all other 
time relationships, are completely identical to those obtained 
above, and, accordingly, the phenomenon of "superresolution" may 
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«Xso Accur here. 


Conclusion 


Proposed in the present study is a new method for processing 
space pulse radio signals for the study of their internal 
structure with a high level of time resolution. The proposed 
method is based on the utilization of the effect of the Fourier 
transform of the envelope of the radio pulses in the interstellar 
plasma, which takes place with the observance of definite 
relationships between the parameters of the pulse and the medium. 

The conditions of applicability of the method are formulated, 
and requirements are also substantiated for the receiving-recording /17 
equipment, which proved considerably lower than for the processing 
methods being utilized at the present time. It is shown that the 
proposed method makes it possible to obtain the autocorrelation 
function of the envelope of the emitted radio pulses with "super- 
resolution”, i.e. , with better time resolution than in the 
recording equipment. 

The obtained results may also find use for problems of the 
synthesis and processing of signals in space communications systems. 

The author would like to express his special gratitude to 
S. V. Pogrebenko for his discussion of an entire series -of questions, 
and also to N. S. Kardashev and M. V. Popov for displaying their 
interest in the study and their valuable comments. 


15 


REFERENCES 


/18 


1. Frenks, D,, Teoriya signalov [Signal Theory], Moscow, "MIR”, 

1978. 

2. Rabiner, L., Gould, B., Teoriya i primsnenie tsifrovoy obrabotki 

signalov [Theory and Application of Digital Probe ssi of Sig- 
nals], Moscow, "MIR", 1978. 

3. Hankins, T. H., Rickett, B. J. , "Pulsar Signal Processing", 

Methods in Computer Physics , Academic Press, 1975, Vol. 14. 

4. Walt, J. R., "Propagation of Pulses in Dispersive Media", Radio 

Sci. 69D/ 11, 1387-1401 (1965). 

5. El-Khamy, S. E., McIntosh, R. E., "OptimupA Transionospheric 

Pulse Transmission", IEEE Trans, of Ant. Propagat. AP-21/3, 
269-273 (1973) . 


